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ABSTRACT 

We present Spitzer IRS spectroscopy of CO2 ice bending mode spectra at 15.2 /^m toward 19 young 
stellar objects with luminosity lower than 1 L Q (3 with luminosity lower than 0.1 L Q ). Ice on dust grain 
surfaces can encode the history of heating because pure CO2 ice forms only at elevated temperature, 
T > 20 K, and thus around protostars of higher luminosity. Current internal luminosities of YSOs 
with L < 1 L Q do not provide the conditions needed to produce pure CO2 ice at radii where typical 
envelopes begin. The presence of detectable amounts of pure CO2 ice would signify a higher past 
luminosity. Many of the spectra require a contribution from a pure, crystalline CO2 component, 
traced by the presence of a characteristic band splitting in the 15.2 /im bending mode. About half 
of the sources (9 out of 19) in the low luminosity sample have evidence for pure C0 2 ice, and six of 
these have significant double-peaked features, which are very strong evidence of pure CO2 ice. The 
presence of the pure C0 2 ice component indicates that the dust temperature, and hence luminosity 
of the central star/accretion disk system, must have been higher in the past. An episodic accretion 
scenario, in which mixed CO-CO2 ice is converted to pure CO2 ice during each high luminosity phase, 
explains the presence of pure CO2 ice, the total amount of CO2 ice, and the observed residual C 18 
gas. 
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1. INTRODUCTION 

While the st andard star formation mod el with constant 
accretion rate (Shu 19771 : iShu et alJll987[ ) predicts proto- 
stars to have a luminosity higher than 1.6 Lq most of the 
time, young stellar objects (YSOs) with luminosity lower 
than 1.6 L w have long been known (IKenvon et al.l 119901 : 
iGreene et all Il994t iHartmann fc Kenvonl 119961 ). The 
Spitzer Legacy Proj ect, From Molecul ar Cores to Planet 
Forming Disks (c2d: lEvans et al.l2003| ) found that 59% of 
the 112 embedded protost ars have luminosity lower than 
1.6 L ()Evans et al.ll2009h . In fact, Spitzer found a sub- 
stantial number of s ources with the internal luminosity 
even below 0.1 L© ()Young et al.l l2004at iDunham et al.l 
2008), which ha ve been called Very Low Luminosity Ob- 
jects (VeLLOs; Idi Francesco et all 12007). The internal 
luminosity (Li n t) of an embedded protostar is the lumi- 
nosity of the central source, excluding luminosity from 
external heating. Recently, an even lower luminosity 
protostar (Li nt < 0.03 L Q ) was disc overed as a binar y 
companion to IRAM04191+1522 IRS (|Chen et al.H2012l ). 
One explanation for such low lumi nosities is that mass 
accretion is not a constant process (IKenvon et ail 119901 : 
lEnoch et al.ll2009bt IDunham et al.ll2010bD . The low lu- 
minosity sources may be going through a low mass ac- 
cretion stage between accretion bursts, explaining their 
currently low luminosities. If mass accretion is episodic, 
sources with low luminosity may have had much higher 
accretion rates in the past, and thus had higher lumi- 
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nosity. Imprints of the high luminosity stage in low lu- 
minosity YSOs would support the idea of episodic mass 
accretion. 

Evidence of past periods of higher luminosity can be 
found in molecular spectra of the ga s and ice phases 
of lo w luminosity s ources (ILed 120071: IVisser fc: Berginl 
2012). In particular. iPontoppidan et al. (2008) analyzed 
the general shape of the 15.2 /im C0 2 ice bending mode 
spectrum toward 50 embedded young stars. The 15.2 
/im C0 2 spectrum can be decomposed into multiple com- 
ponents, including a pure C0 2 ice component. The pure 
CO2 ice can form by two processes. One is CO2 segre- 
gation out of a CO2-H2O mixture. The other is a distil- 
lation process, in which CO evaporates from a CO2-CO 
mixture, leaving pure CO2 behind. The former process 
occurs at a high temperature (50 - 80 K) and the latter 
occurs at a lower temperature (20 - 3 K). Both pure CO 2 
formation processes are irreversible (Hagcn et al. 1983), 
since the bond between pure CO2 ices makes it the most 
stable phase. Once the pure CO2 ice has formed, it will 
not disappear unless it evaporates. Thus the existence of 
pure CO2 ice provides a "chemical memory" of warmer 
conditions in the past. 

The total amount of C0 2 ice, regardless of form, 
is another indicator of the temperature history. A 
number of studies of absorption against background 
stars have show n that CO2 ice must form in the 
prestellar phase (iBergin et all 120051: [Knez et al.1 120051 : 
I Whittet et all 120091 : iChiar et al.1 l20lTh . The shape of 
the 15.2 /im feature to ward background sta rs varies lit- 
tle from cloud to cloud ([Whittet et al.|[2~009f) . with most 
(~ 85 to 90%) of the C0 2 in the C0 2 -H 2 mixture and 
no evidence for a pure C0 2 component, despite varia- 
tions in the ratio of CO to H 2 ice from cloud to cloud 
(| Whittet et all 120091 ) . These ices arise from reactions 
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among H, O, and CO on grain surfaces (e.g.. iChiar et al.l 

USB- 

Observations of the CO2 ice feature toward YSOs 
tend to show higher values of jV(COo ) than the back- 
ground stars ([Pontopp idan et al.ll2008D . even for back- 
gro und stars with rou ghly the same extinction (see Fig. 
5 in lCook et aL|[20Tl . The ratio of C0 2 ice to H 2 ice 
in spectra toward YSOs is also higher and more vari- 
able than toward background stars while the ratio of CO 
to H2O ice is lower for the YSOs, indicative of further 
processing occuri ng during the formation of the YSOs 
(jCook et al.ll27)Tl . 

The total amount of CO2 ice formed during the YSO 
phase depends on the accretion scenario. The episodic 
accreti on scenario gives multip le long periods of low lumi- 
nosity (jDunham et al.ll201QpD instead of the very short 
period of low l uminosity t hat the continuous ac cretion 
model predicts ()Shulll977t lYoung fc Evans! 120051 ). As a 
result, episodic accretion provides more time for CO to 
freeze-out and form CO2 ice. Then, the CO ice can evap- 
orate during episodes of higher accretion and luminosity, 
leaving the C0 2 ice behind. 

The gas phase CO is the final indicator that we con- 
sider. Laboratory experiments show that the binding en- 
ergies of CO and N2 onto dust grain surfaces are similar 
(|Oberg et al.ll200aiBisschop et al.ll2006[ ). After including 
the updated bin ding energy, chemo-dynamical models 
(|Lee et al.ll2004f ) predict more CO, as measured by rare 
isotop es like C 18 Q, th an is observed (e.g., IChen et all 
120091 : iKim et al.l 120111 ). Tying up some of the car- 
bon in CQ 2 ice could help to match the observations 
(|Visser fc Berginl [20121 ). 

To explain the observed gas phase molecular lines and 
the total CO2 ice absorption observed toward one s pecific 
low luminosity source, CB130-1-IRS1. IKim et all poTTf) 
added some sim ple ice reactions to an episo dic accretion 
chemical model (|Lee et al.ll2004t lLeell2007| ). During low 
luminosity periods, the gas phase CO freezes into CO ice, 
and some CO ice turns into CO2 ice, creating a CO-CO2 
ice mixture; then, during high luminosity phases, the CO 
evaporates, leaving behind CO2 ice, which evaporates 
only at a still higher tempera ture (T > 20 K), depend - 
ing on the structure of the ic e (iPontoppidan et al.l i2008'). 
This scenario, de veloped by [Kim et al.l ([201 ID and dis- 
cussed further bv lVisser fc Bergi n (2012), could explain 
the current low luminosity, the total amount of CO2 ice, 
and the strength of the gas phase molecule emission to- 
ward CB130-1-IRS1. 

In this paper, we test the idea furt her with higher res- 
olution Infrared Spectrograph fIRS; iHouck et all I2004J) 
spectra for an enlarged source sample and complemen- 
tary spectroscopy of gas-phase C ls O at 1.3 mm. We ob- 
tained Spitzer IRS, Short-High (SH) mode spectroscopy 
toward 19 young stellar objects (YSOs) with luminosity 
lower than 1 Lq, 4 of them with luminosity lower than 0.1 
Lq. The envelopes around YSOs with luminosity lower 
than 1 L/0 have temperatures lower than the CO evapo- 
ration temperature. The higher resolution spectra allow 
decomposition of the 15.2 /mi CO2 ice bending mode 
spectrum into 5 different component s, including a pure 
CO2 c omponent, using the method of lPon toppida n et al.1 
(2008). In Section [21 we explain the sample selection in 
the current study. We describe the observations in Sec- 



tion [31 In Section [H we describe the analysis of the 
Spitzer IRS spectrum to get the optical depth and the 
column density. In Section[5]we present results of the de- 
composition of ice components. In Section [6l we present 
the setup and results of the chemical models for differ- 
ent accretion scenarios. In Section [71 we summarize our 
findings. 

2. SOURCES 

The sources were selec ted from the lis t of low luminos- 
ity objects compiled by iDunham et al.l ([2008). Table Q] 
includes the source identification numbers and inter- 
nal luminosities from that reference, the Spitzer source 
name, the cores or clouds with which a source is associ- 
ated, Spitzer positions, and the Program ID (PID) of the 
Spitzer IRS observation. The internal luminosities are 
determined by using the correlation betw een 70 /jm flux 
and t he internal luminosity of protostars ([Dunham et al.l 
2008). The sources are Class 0/1 objects with luminosi- 
ties in the ra nge between 0.08 - 0.69 L ^. We also use 
sources from iPontoppidan et al.1 (j2008f ). most of which 
have luminosities between 1 and 10 Lq, but some are 
massive sources with luminosities up to 10 5 L Q . With 
the two samples, we can study objects with luminosities 
from 0.1 L to 10 5 L . 

3. OBSERVATIONS 

The Spitzer IRS observations of 17 sources were ob- 
tained between November 2008 and May 2009 (PID 
50295; PI Michael Dunham) with their background im- 
age. Spectra of two more sources, source number 003 
and 038, were obtained from Spitzer archival data, ob- 
served in 2005, October and November respectively (PID 
20604; PI Adwin Boogert). Background images were not 
observed for these two sources. The observing mode was 
SH, and the resolving power (A/AA) was 600. The data 
were reduced following the IRS pipeline version S18.7 
which produces Basic Calibrated Data (BCD) files. Then 
we used the bad pixel masks to clean the masked pix- 
els in a set of data. We used the Caltcch High-rcs IRS 
pipeline (CHIP) to reduce the 17 observed data with 
background exposure. CHIP is an IDL reduction package 
for reducing high signal t o noise SH Spitzer IRS spectra 
(Pontoppi dan et al.ll2010D . CHIP produces better signal 
to noise level but requires the dedicated background ex- 
posures. For the t wo sources without background images, 
we us ed SMART (jHigdon et alj 120041 : ILebouteiller et all 
2010) to do the background subtraction from their own 
image and to extract spectra from the BCD files. Then 
we trimmed the end of each order of a spectrum. All 
19 sources have clear detections of the CO2 ice bend- 
ing mode absorption spectrum. The overlapping regions 
between orders match well. 

The gas-phase molecular line observations were ob- 
tained at the CSO 6 in 2010 July. The C 18 (J = 2 — >■ 1) 
lines were obtained with the fast fourier transform spec- 
trometer (FFTS) having 8192 channels with a 500 MHz 
total bandwidth. We used the 230 GHz heterodyne re- 
ceiver in position switching mode. The main beam effi- 
ciency was 0.8, and the velocity resolution was 0.17 km/s. 

6 The Caltech Submillimeter Observatory was supported by the 
NSF. 
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4. ANALYSIS 

4.1. Reduction of the IRS Spectra 

The 15.2 /mi CO2 bending mode spectrum is located 
between the broad 9.7 /mi, and 18 - 20 /mi silicate fea- 
tures, and near the 12 /mi H 2 band. To remove the 
baseline, we constructed a third-order polynomial in the 
range of 13.3 — 14.7 /im and 16.2 — 19.3 /im,. and added 
a gaussian feature centered at 608 cm -1 with FWHM of 
73 cm" 1 for the blue wing of th e 18 - 20 /im silicate 
feature (|Pontoppidan et al .1 l2008h . The 15.2 /mi C0 2 
bending mode spectra are plotted in Fig.Q]- Fig. [5] The 
best fit baselines are given as red dashed lines in the left 
panels of Fig. Q] - Fig. [5] After the continuum subtrac- 
tion, we calculated the optical depth of each source by 

r = — In (^j^j > where Iq is the continuum intensity. Op- 
tical depths are plotted in the right panels of Fig. Q] - 
Fig. G3 

4.2. Decomposition into Ice Components 

We used a number of laboratory spectra of CO2 
ices to analyze the 15.2 /mi CO2 bending mode. The 
available laborator y spectra are CCyB^O mixtures, 
CQ2: CO mixtures (lEhrenfreund et al.l ll997l ). and pure 
CO2 (|van Broekhuizen et al.l 120061 ). A minimum of five 
unique components has been requir ed to fit all previ- 
ously observed 15.2 /im C0 2 spectra (jPontoppidan et aTl 
\mm . They are a dilute C0 2 (C0 2 :CO = 4:100), an 
H 2 0-rich component (CC^H^O = 14:100), a mixture of 
CO and C0 2 (C0 2 :CO = 26:100 and C0 2 :CO =70:100), 
a pure CO2 component, and a shoulder component. We 
did not include the shoulder component, since the shoul- 
der component appears only in massive young stars, but 
we use both C0 2 :CO = 26:100 and C0 2 :CO =70:100 
mixtures as separate components. Thus, the number 
of CO2 ice components is still 5 in the current study. 
We use the laboratory data with the grain shape of 
the continuous distribution of ellipsoids (CDE) to con- 
vert optical constants to opacities. This is a method 
to model irregularly shape d grains, successfully used 
in IPontoppidan et al.l (|2008l) . The laboratory ices have 
been obtained at temperatures between 10 K and 130 K. 
The laboratory data are not strongly dependent on the 
temperature. For example, the strength of the pure CO2 
ice component decreases by 1% when the temperature in- 
creases from 10 K to 50 K. When temperature increases 
to 80 K, the strength of the component decreases by 9%. 
The envelopes of low luminosity sources have tempera- 
tures lower than 20 K. We use the laboratory ices with 
a temperature at 10 K, except that the pure CO2 com- 
ponent is at 15 K, for which the laboratory data are 
available. 

The laboratory data for the H^O-rich CO2 ice compo- 
nent lies on top of a broad 12 /im H 2 ice feature. So 
we fitted a baseline with a third order polynomial, and 
then subtracted the baseline from the laboratory data. 
The resulting laboratory data used to decompose CO2 
ice features are plotted in Fig. [6] The dilute CO2 com- 
ponent is located at the red part of 15.2 /mi bending 
mode. The C02:CO = 70:100 component is located at 
the blue part of 15.2 /mi bending mode. Only the pure 
CO2 ice component shows a double peaked feature. 

We varied only the strength of each component to fit 



the observed optical depth of the CO2 feature. Then we 
used a x 2 minimization technique to obtain the best fit. 
We have 5 free parameters, the strengths for each of the 
5 different ice components. The reduced \ 2 of the best fit 
is listed in Table [2j The best fit laboratory components 
are plotted separately and in total in the right panels 
of Fig. [T]- Fig. [5] In all sources, the water-rich CO2 ice 
dominates the shape and the strength of the spectra. The 
red wing part of the spectrum is determined by the water- 
rich component. The peak strength is determined by the 
sum of all the components. The double peaked features, 
unique to the pure CO2 ice, generally agree in wavelength 
with the laboratory data. We did the component analysis 
with and without the pure CO2 ice component included 
in the fit. 

4.3. Column Densities of Ices 

To obtain a column density from the optical depth, we 
integrated the optical depth over wave number, and then 
divided it by the band strength, A: 



N(cm 



r(u)dv/A. 



(1) 



We used the band strength from lGerakines et al.l (|1995D . 
For the pure CO2 component, the C02:CO = 4:100 com- 
ponent, and the C02:CO = 26:100 component, A — 
1.1 x 10" 17 cm molecule" 1 . For C0 2 :CO =70:100, A = 
0.97 x 10~ 17 cm molecule" 1 , and the component with 
C0 2 :H 2 = 14:100 has A = 1.57 x 10" 17 cm molecule" 1 . 
The column densities for each component are listed in 
Table H 

4.4. Column Density of Gas Phase CO 

The CSO data were reduced with CLASS 7 . The molec- 
ular line results are summarized in Table [3J The column 
density of C 18 is obtained from 



N. 



c 18 o 



3kQ exp(Ej/kT ex ) 
8ir 3 vfi 2 J 



Tr dv, 



(2) 



where Ej = hBJ(J+l), B is the rotational constant, T ex 
is the excitation temperature, /i is the dipole moment, 
and Q is the partition function. The parameters B and 
fi for C 18 come from Table 3 of lLee et al.l (|200l . We 
assumed T ex = 10 K to derive the column density. 

5. RESULTS 

Double peaked features are clearly seen for six sources 
(30%), as indicated in Table [2j As Fig. [6] indicates, only 
the pure CO2 ice component has a double peak at wave- 
lengths 15.14 /tm and 15.26 /mi. All the sources with 
obvious double-peaked features have fits for column den- 
sity of the pure CO2 component that are more than 10 
tr, and the value of x 2 improved by more than 20% when 
the pure C0 2 ice component was included. There are 
three other sources without obvious double-peaked fea- 
tures that satisfy the same criteria for signal-to-noise and 
improvement in x 2 , also marked in Table [5J We also con- 
sider these to be detections, for a total of 9 detections out 
of 19 sources (47%). In fact, all but 3 sources have best 
fit values of CO2 ice column densities with signal-to-noise 

7 http: / /www.iram.fr /IRAMFR / GILDAS 
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greater than 10, but we require the improvement in x 2 
because uncertainties in a multi-component fit may be 
unreliable. 

We detected the pure CO2 ice component with a signif- 
icant double peak even from t he source with L int = 0.08 
L© (Source number: 003). iPontoppidan et all (|2008[ ) 
found the pure CO2 ice in 40% of the total 56 sources 
with clear detections in sources with luminosity higher 
than 1 L©. The current result shows that the presence 
of pure CO2 ice is at least as common in low luminosity 
protostars as it is in high luminosity protostars. This 
is all the more remarkable because absorption in the 
surrounding molecul ar cloud from unprocessed CO2 ice 
(jWhittet et al.ll2009D will make the pure CO2 ice signal 
harder to see. 

CO evaporates from CO2-CO ice mixture at a temper- 
ature 20 - 30 K. We take the minimum required tem- 
perature for pure CO2 ice formation as 20 K. It is be- 
yond the scope of this paper to explore the parameter 
space of the structures of the protostellar envelopes, but 
through model studies, the envelope inner radii of em- 
bedded protostars are > 100 AU. (Serp ens 1: 600 AU 
Enoch et al.ll2009aL IR AM 04191: 140 A UlD unham et all 
20061 CB130-1 : 350 A UI Kim et al.ll20lH L673-7: 140 AU 
Dunham et al.ll2010aft . Many of these individual studies 
used simple 1-D or 2-D geometry to approximate the 
complicated transition between the envelope and disk, 
and the envelope inner radius is a way of parameterizing 
this transition with a single value. To match observed 
SEDs, one must have an inner cut-off to a spherical den- 
sity profile to avoid having extremely high density mate- 
rial, and thus optical depth. Though 100 AU is not an 
exact value representing the inner radii of all envelopes, 
it is a minimum value compared to values found in many 
models. 

Using DUSTY (|Ivezic et al.l Il999f ). we calculated a 
dust temperature around sources with luminosities of 
0.7 Lq (the highest current luminosity in our sam- 
ple) and 150 L^, more c haracteristic of the sample of 
Pontoppi dan et al.l (|2008t ). For these models, we de- 
creased the inner radius so that we could trace the dust 
temperature to smaller radii. The dust temperature for 
the L int = 0.7 L Q protostar (Fig. [7]) is lower than 20 K ev- 
erywhere outside about 400 AU. A detailed calculation of 
the chemistry preceding and during the Firs t Hydrostatic 
Core (FHSC) stage confirms this result (jFuruva et al.l 
2012). In contrast, the model with 150 L© has an ex- 
tended region (r < 8000 AU) with dust temperature 
above 20 K. If the luminosities had always been less than 
or equal to the current values, pure CO2 ices should be 
prominent only in the high luminosity sample, contrary 
to our observations. 

We can also make a quantitative comparison of the 
curren t sample and the sample from IPontoppidan et al.l 
(2008). We plot column densities of each CO2 ice com- 
ponent versus source luminosities in Fig. [5] We divide 
the sources for convenience into three groups: the cur- 
rent, low-luminosity [Lint < 0-7) sample; the sources 
from IPontoppidan et al.l (|2008D with intermediate lumi- 
nosity (0.7 < L int < 100 L Q ); and the high-luminosity 
(L int > 10 3 L Q ) sources. The sources in the current 
study, with luminosity lower than 0.7 L Q , have a simi- 
lar pure CO2 column density range as the intermediate- 



luminosity sources. For example, the maximum pure 
CO2 ice column density detected toward sources with 
luminosity higher than 1.0 L Q is 1.06 x 10 18 cm~ 2 , while 
the maximum pure CO2 ice column density is 6.8 x 10 17 
cm~ 2 in the source with internal luminosity 0.54 Lq. The 
amount of pure CO2 ice does not depend on the current 
luminosity. 

Table [land Fig. [8] show that the H 2 0-rich C0 2 ice 
component is dominant in all sources. Most of the CO2 
ice exists as a mixture with water regardless of the source 
luminosity. It takes more energy to segregate CO2 ice 
and H2O ice, so H^O-rich CO2 ice is more stable than 
CO-CO2 ice mixtures. In the low-luminosity sources, the 
mixture with CO exhibits more scatter and higher val- 
ues in some cases than in the high-luminosity sources. 
The sources with high luminosities may have evaporated 
much of the CO from CO-CO2 ice leaving pure CO2 
behind. The total CO2 ice amount is also low in the 
high-luminosity sources compared the the spread seen in 
the low-luminosity sources. The high-luminosity sources 
could evaporate all the local CO2 ices, leaving only the 
ices in the extended cloud. The pure C02/total CO2 ra- 
tios have almost the same range of variation in the low 
and intermediate-luminosity groups. This result adds to 
the evidence for higher luminosities in the past for the 
sources with current low luminosities. 

6. DISCUSSION 
6.1. Chemical Models with Episodic Accretion 

One plausible reason for the presence of the pure CO2 
ice compo nent in low luminosity young stars is episodi c 
accretion (jKenvon et al.l 119901: iVorobvov fc Basul 120051) . 
In this model, the material falling in from the envelope 
is not continuously accreted onto the star. Instead, it 
piles up on a disk and accretes onto a star in a short- 
lived burst, when the disk becomes unstable. The burst 
period can give enough heating to form pure CO2 ice in 
the envelope. Once the pure CO2 ice has formed, it will 
persist. A series of "freeze-thaw" cycles combined with 
irreversible chemical processes can make distinct chem- 
ical signatures not possible in a model with monotonic 
temperature behavior. 

In addition to the pure CO2 evidence, the total amount 
of CO2 ice, including that in mixed ices, can be explained 
by long periods of low luminosity between episodic ac- 
cretion bursts, as predicted in an episodic accretion sce- 
nario. The source luminosity steadily increases in a con- 
tinuous accretion model as the mass of the central star 
grows. Larger and larger fractions of the envelope be- 
come too warm for freeze-out. In contrast, the episodic 
accreti on model has multiple long periods of low lumi- 
nosity (jDunham et al.l I2010bf) , punctuated by episodes 
of higher luminosity. More ice can form during the low 
luminosity periods. 

With chemical evolution modeling, we can relate the 
amount of observed CO2 ice to the accretion scenario. 
We use the evo lutionary chemo-dynamical model by 
ILee et all (|2004f) to test this idea. The model calcu- 
lates the chemical evolution of a core from the prestellar 
core to the embedded protostellar core stage. At each 
time step, the density profile, the dust temperature, the 
gas temperature, and the abundances are calculated self- 
consistently. 
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The original chemical models can produce CO2 ice 
only by freezing of gas-phase CO2. We added a 
second route to CO2 ice production into the chemi- 
cal network. Until recently there was no known ex- 
act chemical p athway for CO2 ice form a tion in cold 
environments {D'Hendecourt et aD 119851: iRoser et aL 
[200llRuffle fc Herbstll2001fl . Recently. IGarrod fc Raurv 
( 20111) suggested that CO ice turns into CO2 ice on 
cold interstellar dust grain surfaces in the presence of 
OH radicals. Lab oratory experiments support this idea 
(|Oba et al.l l2010f ) . The process can occur at dust tem- 
peratures as low as ~ 10-12 K. In our study, we did not 
include the full reaction of CO and OH radicals because 
the kinetics of CO2 ice formation on surfaces is an active 
field of research with no clear answer so far. Instead we 
adopted a simple model. When the dust temperature is 
below the CO ice freezing point, (100 — x)% of the CO 
gas turns into CO ice, and x% turns into CO2 ice. 

We ran two kinds of chemical models, those with and 
those without the additional pathway to CO2 ice from 
CO. For each kind of chemical model, we ran two kinds 
of models for the luminosity evolution, for a total of four 
kinds of models. The first luminosity model assumes 
continuous accret ion. We adopt t he lum inosity evolu- 
tion profile from lYoung fc Evans! (120051) including the 
First Hydrostatic Core stage (lLarson"fl9 69: Boss 119931 : 
iMasunaea et alj|1998t iFuruva et al.ll2012D . Without the 
FHSC, the luminosity evolution does not have a stage 
with luminosity lower than 1.0 Lq, so we use the model 
with the FHSC stage. In this picture, the low-luminosity 
sources would be in the FHS C stage. The second i s an 
episodic accretion model from lDunham et"ail (|2010b[) . In 
this model, the episodic accretion is included in a simple, 
idealized way. The model assumes no accretion from disk 
to star most of the time, but steady infall from envelope 
to disk. When the disk mass reaches 0.2 times the stel- 
lar mass, then an accretion burst occurs with luminosity 
100 times more than the u sual value. We used the model 
of a 1 M envelope from iDunham et al.l (|2010bD . The 
luminosity evolutions of the continuous model and the 
episodic accretion model are plotted in the top panels of 
Figs. [9] and [TUl respectively. At early time steps, up to 
20000 yrs, both the continuous accretion model and the 
episodic accretion model evolve in the same way, going 
through a FHSC (First Hydrostatic Core) stage. The 
luminosity increases monotonically until it reaches 0.01 
Lq. After the FHSC stage ends, the accretion luminos- 
ity increases rapidly, so there is a luminosity gap between 
0.01 L Q and 1 L Q in the continuous accretion model and 
a gap between 0.01 L© and 0.1 L in the episodic accre- 
tion model. 

At each time step, we calculat e the dust temperature 
with DUSTY (llvezic et aflll999D and the gas tempera- 
ture with a gas en ergetics code (|Dotv fc Neufeldl fl997t 
Young et al1l2004b[) using the same modeling set up as 
Kim et al. Then the chemical evolution of 512 



gas parcels is calculated as each falls into the central re- 
gion. We assumed the surface binding energy of species 
onto bare Si02 dust grains. The chemical calculation in- 
cludes interactions between gas and dust and gas-phase 
reactions. Then the abundance profile at each time step 
is calculated. To compare models with ice observations, 
we calculated the total CO2 ice column density at each 
time step. We do not track the different components of 



the CO2 ice, as the chemical reactions between them are 
not well established, so we will compare the results to 
only the total CO2 ice. 

To compare to the observations of gas phase CO, we 
use the distributions of density, gas temperature, abun- 
dance, and velocity at each time step to calculate the 
distribution of C ls O over its energy levels and then sim- 
ulate the observations by integrating through the enve- 
lope and convolving the intensity with a beam matched 
to the observations. Then we derive the "observed" col- 
umn density of C 18 0, using the same procedure as used 
for the observations. The column density of gas phase 
CO is then obtained by multipl ying the C 18 Q colum n 
density by the isotope ratio (540: |j0rgensen et al.ll2004j) . 
Because of the limited spatial resolution of the obser- 
vations and assumptions involved in calculating the ob- 
served column density, this "observed" column density 
from the models will differ from the true column density, 
but it provides the best comparison to observations. 

Models were run with different values of x, the per- 
centage of freezing CO that goes into CO2 ice. The best 
match to the observations was found for x = 10, so we 
focus the rest of the discussion on those models. First 
we discuss the model with a constant accretion rate, fol- 
lowed by that with episodic accretion. 

In the model with constant accretion (Fig. [SJ the lu- 
minosity increases monotonically with time once the first 
hydrostatic core forms (top panel). The other two pan- 
els in Figure [9] show the total column density in CO2 
ice and in CO gas; they include the prestellar phase 
(t < 0), when the core is assumed to pass through a series 
of Bonnor-Ebert spheres of increasing central concentra- 
tion. Since evolution in this phase is slow, the time axis is 
condensed. During this phase, the "observed" gas-phase 
CO increases by a factor of about 2 because the increas- 
ing central density increases the excitation of the C 18 0. 
The CO2 ice column density increases as gas-phase CO2 
freezes out. These processes accelerate at the later stages 
of highest central density, but the freeze-out of CO be- 
gins to decrease the gas phase column density toward the 
end of the prestellar evolution. At t = 0, the prestellar 
core evolution is ended by the formation of the FHSC. 
At t > 0, the luminosity of the FHSC grows slowly, and 
some CO2 ice evaporates, but it still dominates the gas 
phase CO column density until the true protostar forms 
at t = 2 x 10 4 yr. There is a slight increase in the CO2 
ice toward the end of the FHSC phase owing to the con- 
version of CO ice to CO2 ice. After t — 2 x 10 4 yr, the 
luminosity increases as the stellar mass grows (L oc M) 
and the CO2 gradually declines while the gas phase CO 
increases. However, the CO2 ice column density is still 
about twice the gas-phase CO column density. 

Fig.fTOlshows the same graphs for the model of episodic 
accretion. The evolution through the FHSC phase is 
identical to that of the constant accretion model, but 
the differences begin with the first episode of enhanced 
accretion. The dramatic increase in luminosity causes 
a sharp increase in the gas phase CO, at the expense 
of the CO2 ice. After the accretion burst ends, the CO 
begins to freeze again and the CO2 ice builds up. This 
cycle is repeated, with the balance of the C shifting be- 
tween ice and gas through each accretion cycle, but more 
ice forms in later quiescent cycles because there is more 
time between bursts and the luminosity between bursts is 
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slightly lower as the steady state accretion rate is lower. 
Also, some CO2 ice persists in the cooler parts of the en- 
velope in each cycle. After the last accretion burst, the 
luminosity declines slowly and the CO2 ice increases at 
the expense of the gas phase CO, but this process is lim- 
ited by processes that continually evaporate some CO. 
The evolution o f the ice and gas comp onents is similar 
to that shown in lVisser fe Berginl (|2012l ); differences can 
mostly be attributed to the fact that our model follows 
the gas as it fal l s in, unlike the single-point models in 
iVisser fe Berginl (pOTl ). 

6.2. Comparison to Observations of Total CO2 Ice 

In Fig. QTJ we compare the model results and the ob- 
servational data. The green dots are the observed total 
CO2 ice column density fr om both the current study and 
iPontoppidan et ail (|2008l ). The observed source lumi- 
nosities range from 0.1 L© to 90 L©. The continuous 
model does not have luminosities higher than 2.4 L©, 
and the simple, bimodal, episodic accretion model does 
not cover the luminosity range between 2.6 L© and 90 
L©. Therefore, we compare the model and observations 
only between 0.1 L© and 2.5 L©. 

Both the original episodic accretion model and the con- 
tinuous accretion model, with the only path to CO2 ice 
being freeze-out of CO2 gas, underestimate the observed 
column densities by an order of magnitude. The detailed 
mode l of FHSC formation and evolution (jFuruva et al.l 
120121) confirms this result. Matching the observations 
requires the extra pathway in which CO ice is converted 
to CO2 ice. Of the two models with this pathway, the 
continuous model predicts only a small range of CO2 
ice column densities. The episodic accretion model with 
conversion of CO ice to CO2 ice describes both the lu- 
minosity spread and the spread in observed total CO2 
column density reasonably well. 

The total column density of CO2 ice at the end of the 
prestellar core evolution is about 1 x 10 18 cm" 2 , is similar 
to column densities seen toward background stars (e.g. 
IWhittet et al.1120091) . Consequently, the accumulation of 
CO2 ice in the quiescent cloud, before the formation of 
the core would roughly double the column density in our 
models, which do not explicitly include this evolution. 
Differences from cloud to cloud and from YSO to YSO 
in the amount of this absorption from the quiescent cloud 
could explain some of the scatter, allowing the continuous 
accretion model to fit the data better. However, much of 
the scatter lies in the regime of iV(C02) above the val- 
ues u sually seen toward background stars ()Whittet et al.l 
[20091) . 

6.3. Comparison to Observations of C 18 O Gas 

Another check on the models is provided by the resid- 
ual amount of C 18 gas. Leaving out the two sources 
with the highest CO column density, the average C 18 
gas column density is 9.18 x 10 14 cm" 2 . We round off 
to ~ 10 15 cm" 2 as a residual C 18 gas column density, 
which is plotted as a horizontal line in Fig. [12l The 
four models are plotted in four different panels for com- 
parison to the data. After a certain amount of time, 
the C 18 gas predicted by the constant accretion mod- 
els converges to a roughly constant value, N ~ 6 x 10 14 
cm" 2 , lower than the mean value. In the episodic ac- 
cretion model, C 18 gas evaporates after each episodic 



accretion burst, but N ~ 10 15 cm -2 most of the time, as 
observed. Of course, C 18 emission can also come from 
the larger cloud. An example may be source 182, with 
the largest C 18 column density; this source is in the 
Ophiuchus cloud, which has large column densities even 
along lines of sight that do not intercept YSOs. 

The different predictions arise from the different abun- 
dance profiles. The abundance profiles at time step 60000 
yr, the last time step in the continuous accretion model, 
are plotted in Fig. [T3J along with those from the episodic 
accretion models, which are in a high luminosity state 
at this time. The high luminosity has heated the enve- 
lope, and thus raised the gas phase CO abundance out to 
much larger radii, compared to the continuous accretion 
model. 

The column density versus luminosity graphs (Fig. I14[) 
compare the models to observations. Seven sources 
are scattered around the predictions of the model with 
episodic accretion and CO to CO2 ice conversion. Two 
sources lie well above and two lie well below the predic- 
tions. The source with the least amount of C 18 gas is 
CB 130-1. That source was best fit with a model in whic h 
80% of the CO ice turns into C0 2 ice (|Kim et al.ll2011l ). 
which is much higher than used here. Most source are 
better matched if 10% of CO ice turns to CO2 ice. 

The current models provide a proof of concept for 
episodic accretion models to explain the observed CO2 
ice and gas-phase CO abundances, but they can be im- 
proved by future work. The episodic accretion included 
in this study is a simple and idealized model with ex- 
treme luminosity variations and does not include the ab- 
sorption from the larger cloud. It does not cover the 
luminosity range between 2.3 L© and 90 L©. Fig. [TT] 
shows that there is a luminosity gap between 2.3 L© 
and 90 L© in the episodic accretion scenario, while the 
observed sources do have luminosities in the range be- 
tween 3 L© and 100 L©. Using more realistic models 
of the accretion history and a range of envelope m asses 
(jOffner fe McKedl2T)dll iDunham fc Vorobvovll20T2l ) will 
allow more detailed comparison to the observations. 

6.4. Is Episodic Accretion a Unique Explanation? 

In this paper we have s hown that a scenario of episodic 
accretion as modeled by IDunham et al.l (|2010bl ) can ex- 
plain the presence of a component of pure CO2 ice in 
the observed 15.2 /im absorption feature, the total col- 
umn density of all CO2 ice, and the column density 
of gas-phase C 18 in a sample of low-luminosity pro- 
tostars. Coupled with the ability of such a scenario to 
also exp lain the observed protos t ellar luminosity distri- 
bution dVorobvov fc Basul [20091: IDunham etahl l2010bl: 
IDunham fc Vorobvovl 120121 ) and with the existence of 
other, indirect observational evidence for variable accre- 
tion (Dunham & Vorobyov 2012 and references therein), 
there seems to be strong evidence in favor of the existence 
of episodic accretion. 

However, models predicting other mass accretion and 
luminosit y histories have also been proposed. As first 
noted bv iKenvon et al.l (jl990l ). a model in which accre- 
tion rates peak early and then decline through the re- 
mainder of the protostellar phase could potentially match 
the observed protostellar luminosity distribution. Mod- 
els of this ty pe have been advanced by a number of 
authors (e.g.. kndre et al.l 119991: lOffner fc McKed 120111: 
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Hansen et all |2012T >. iDunham fc Vorobvovl (|2012T > re- 
cently argued that such a scenario by itself is insuf- 
ficient to match observed protostellar luminosities, al- 
though their conclusions were based only on the specific 
set of models considered in their study. A model with 
an initial peak in the luminosity followed by a long pe- 
riod of decline might reproduce the observations if the 
initial peak were high enough to produce sufficient pure 
CO2 ice and the time spent at lower luminosities were 
long enough to convert sufficient CO ice to CO2 ice. 
Full chemical modeling of such an accretion scenario is 
required to fully evaluate how well it could reproduce 
our observations, but this is beyond the scope of our 
present study. While we acknowledge that such a sce- 
nario could possibly explain our ice and gas-phase ob- 
servations, we consider it somewhat unlikely given the 
existence of sources like IRAM04191+1522 (source 003) 
that clearly show double-peaked features indicative of 
pure CO2 ice. These are still deeply embedded, Class 
sources, too young to have already spent sufficient time 
in the low-luminosity, dec lining accretion ph ase in typical 
models of this type (e.g.. IAndre et al1ll999t ). 

Other accretion scenari os may also expla i n obs erved 
protostellar luminosities. Offner fc McKed (20 pre- 
sented analytic derivations of the protostellar luminos- 
ity distribution for different models and concluded that 
models that tend toward a constant accretion time rather 
than constant accretion rate produce a greater spread 
in luminosi ties and are i n bett er agreement with ob- 
servations. lHansen et al.l (|2012f ) presented simulations 
including protostellar feedback, both protostellar radia- 
tion and outflows, and argued that such simulations can 
match the means, medians, and standard deviations of 
the observed protostellar luminosity distribution with- 
out episodic accretion. However, they did not compare 
the full observed and modeled lumin o sity d istributions. 
The accretion rates in lHansen et al.l (120121) and a sub- 
set of the models in lOffner fc McKed (|2011l ). the subset 
they called the "tapered" models, feature an early peak, 
followed by a slow decline in the accretion rate through- 
out the remainder of the embedded phase. As discussed 
above, such models could possibly reproduce our obser- 
vations, although such models will also be challenged to 
predict the existence of very young, deeply embedded, 
very low luminosity protostars with significant pure CO2 
ice. 

Detailed calculations of the chemistry for specific evo- 
lutionary histories will be needed to test these other mod- 
els. One possible problem with models that feature one 
period of high luminosity followed by a long period of low 



Andre, P., Motte, F., & Bacmann, A. 1999, ApJ, 513, L57 
Bergin, E. A., Melnick, G. J., Gerakines, P. A., Neufeld, D. A., & 

Whittet, D. C. B. 2005, ApJ, 627, L33 
Bisschop, S. E., Fraser, H. J., Oberg, K. I., van Dishoeck, E. F., 

&; Schlemmer, S. 2006, A&A, 449, 1297 
Boss, A. P. 1993, ApJ, 410, 157 

Chen, J.-H., Evans, N. J., Lee, J.-E., & Bourke, T. L. 2009, ApJ, 
705, 1160 

Chen, X., Arce, H. G., Dunham, M. M., & Zhang, Q. 2012, ApJ, 
747, L43 

Chiar, J. E., Pendleton, Y. J., Allamandola, L. J., et al. 2011, 
ApJ, 731, 9 



luminosity is that CO ice will form over the pure CO2 
ice, rendering it less visible. 

Ultimately, our results do not prove that episodic ac- 
cretion occurs. They do, however, demonstrate that the 
past accretion rates must have been higher than the cur- 
rent rates for about half of our targets, and they do show 
that a scenario of episodic accretion is fully consistent 
with our observations. 

7. CONCLUSION 

We detected the pure CO2 ice double-peaked feature 
from 6 out of 19 low luminosity sources with L < 1 L Q . 
Fits to the absorption profile provide strong evidence for 
pure CO2 ice in 3 more sources, for a total of 9 (47%). 
The minimum required temperature to form pure CO2 
ice is 20 K, which is not found in any substantial part of 
the envelope at the current luminosities of these sources. 
Detection of pure CO2 ice strongly indicates higher lumi- 
nosities in the past. During the time of higher luminosity, 
the pure CO2 ice forms and persists into the current low 
luminosity stage. The existence of pure CO2 ice toward 
such low-luminosity sources is very strong evidence of a 
phase of higher luminosity in the past, consistent with 
models of episodic accretion. It also presents a challenge 
to any alternative model for explaining the luminosity 
spread seen toward YSOs. 

In addition to th e luminosity spread of sources 
(Dunham et al. 2010b, Dunham and Vorobyov 2011), the 
total CO2 ice amount and the residual gas phase CO are 
best explained by episodic accretion with conversion of 
CO to C0 2 ice. If 10% of CO goes into C0 2 ice when the 
gas phase CO freezes, observations of most sources are 
reasonably well matched. There is a substantial spread 
of values of both total CO2 ice and gas phase C 18 in 
the observations, which can be explained in an episodic 
chemistry. However, variations in the ice formation in 
the larger cloud can also contribute to this scatter. De- 
terminations of other ice components and extinction in 
the surrounding cloud can help to distinguish these ef- 
fects. 
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TABLE 1 

The low luminosity sources observed with Spitzer IRS 



Source Number a 


Spitzer Source Name 


c2d Core Region 


RA 


DEC 


Internal luminosity 13 


PID 








(J2000) 


(J2000) 


(L ) 




003 


SSTc2d J042200.41+153021.2 


IRAM04191+1522 


04:22:00.41 


+15:30:21.2 


0.08 


20604 


005 


SSTc2d J044112. 65+254635.4 


TMC1 


04:41:12.65 


+25:46:35.4 


0.36 


50295 


016 


SSTc2d J124539. 96-552522. 4 


DC302.1+7.4 


12:45:39.96 


-55:25:22.4 


0.42 


50295 


017 


SSTc2d J130736.89-770009.7 


DC303.8-14.2 


13:07:36.89 


-77:00:09.7 


0.46 


50295 


021 


SSTc2d J165719. 63-160923. 4 


CB68 


16:57:19.63 


-16:09:23.4 


0.54 


50295 


023 


SSTc2d J171122. 18-272602.0 


B59 


17:11:22.18 


-27:26:02.0 


0.42 


50295 


024 


SSTc2d J181616. 39-023237.7 


CB130-1 


18:16:16.39 


-02:32:37.7 


0.15 


50295 


038 


SSTc2d J212407.58+495908.9 


L1014 


21:24:07.58 


+49:59:08.9 


0.09 


20604 


043 


SSTc2d J222959. 52+751403.1 


L1251 


22:29:59.52 


+75:14:03.1 


0.21 


50295 


060 


SSTc2d J032637.46+301528.1 


Perseus 


03:26:37.46 


+30:15:28.1 


0.69 


50295 


063 


SSTc2d J032738. 26+301358.8 


Perseus 


03:27:38.26 


+30:13:58.8 


0.20 


50295 


068 


SSTc2d J032845. 29+310542.0 


Perseus 


03:28:45.29 


+31:05:42.0 


0.26 


50295 


078 


SSTc2d J032923. 47+313329.5 


Perseus 


03:29:23.47 


+31:33:29.5 


0.20 


50295 


090 


SSTc2d J033229. 18+310240.9 


Perseus 


03:32:29.18 


+31:02:40.9 


0.20 


50295 


092 


SSTc2d J033314.38+310710.9 


Perseus 


03:33:14.38 


+31:07:10.9 


0.14 


50295 


109 


SSTc2d J034421. 36+315932.6 


Perseus 


03:44:21.36 


+31:59:32.6 


0.11 


50295 


124 


SSTc2d J125342.86-771511.5 


Chamaeleon II 


12:53:42.86 


-77:15:11.5 


0.14 


50295 


161 


SSTc2d J160115. 55-415235. 4 


Lupus 


16:01:15.55 


-41:52:35.4 


0.08 


50295 


182 


SSTc2d J162705. 23-243629. 5 


Ophiuchius 


16:27:05.23 


-24:36:29.5 


0.15 


50295 



a The source number comes from lDunham et al.l I200gf) . 

b The source luminosity comes from the relation between 70 flux and internal source luminosity (Dunham ct al. 2008). 



TABLE 2 

Ice column densities of the CO2 component 



Source Number Pure CO2 CO:CO 2 =100:4 CO:CO 2 =100:26 CO:CO 2 =100:70 H 2 O:CO 2 =100:14 Total C0 2 \ 2 



003 a 


1 


.717 


± 





.130 





.063 


± 


0.005 





.305 


± 


0.023 


1 


.336 


± 


0. 


.106 


5 


.949 


± 


0. 


.417 


9 


.370 


± 


0. 


.680 





.034 


005 





.297 


± 





.014 





.000 


± 


0.001 





.801 


± 
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Note. — All column densities arc in 10 17 cm -2 . 
a Clearly double-peaked feature. 

High signal-to-noise column density of pure CO2 feature and \ 2 improves by more than 20% when pure CO2 is included. 
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TABLE 3 

The C 18 ( J = 2 -> 1) observation of low luminosity embedded sources 



g Number 


rms 


f T*. rlv 


VLS ft 


Av 


T*. 


Column Density 




(IS.) 


^iv Kin b j 


(km s ) 


(km s ) 




(cm -2 ) 


003 




1.68 


6.38 


0.76 


2.06 


1.31c+15 


005 


0.138 


1.31 


5.24 


0.70 


1.77 


1.02e+15 


021 


0.144 


1.25 


5.03 
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9.69e+14 
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3.13 


3.45 


1.03 


2.86 


2.44e+15 
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0.48 


7.61 


0.71 


0.64 


3.74e+14 
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0.216 


1.20 


-4.22 


0.82 


1.36 


9.31c+14 
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0.58 
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0.52 
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4.52e+14 


063 


0.135 


1.22 


5.20 


2.59 


0.44 


9.49e+14 
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0.141 


0.99 


8.34 


1.41 


0.66 


7.72e+14 
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0.138 


1.71 


8.96 


0.84 


1.90 


1.33e+15 


161 


0.172 


1.39 


4.02 


0.77 


1.70 


1.08e+15 


182 


0.199 


4.69 


4.72 


1.20 


3.66 


3.65e+15 



Reference 



(1) 
(2) 



Note. — Reference: (1) Chen et al (2012, in prep), f2) IKim et al.l 1 1201 II) 
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Fig. 1. — Left panels : observed flux (black solid line) and best fit continuum (red dashed line). Right panels : optical depth (black 
solid line), sum of all the ice components (red solid line, best fit model), pure CO2 (yellow dash-dot), H2 0-rich CO2 ice component (blue 
dashed), CO-CO2 mixtures (purple dashed). 
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Fig. 2. — Continued 




Fig. 3. — Continued 
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Fig. 4. — Continued 
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Fig. 5. — Continued 
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Fig. 6. — A plot of laboratory data. The yellow line is a pure CO2 ice, the light blue line is a water-rich CO2 ice, the blue line is a 
dilute CO2 ice, and purple lines are CO:C02 ice mixtures. Only the pure CO2 ice component shows a double peak feature. A baseline is 
subtracted from a water-rich CO2 ice component. 
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Fig. 7. — Dust temperatures as a function of radius for two different luminosities, characteristic of the current sample (L = 0.7 Lq) and 
the previous sample (L = 150 Lq). The inner radius in both models is 10 AU, the same as was used in the chemical models used in this 
paper. The envelope inner radii of embedded protostars are about 100 AU. The horizontal dashed line indicates a dust temperature of 20 
K, the minimum required for producing pure CO2 ice. 
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Fig. 8. — The column densities of the different CO2 ice components and t he total plotted versus source luminosities. The red points are 
from the current study; the black points are from Pontoppidan ct al. ( 2008) 
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Fig. 9. — Top panel: The luminosity evolution of the continuous accretion model. Middle panel: The CO2 ice column density evolution 
in the constant accretion model with 10% CO ice to CO2 ice conversion. Bottom panel: The gas phase CO column density evolution in the 
constant accretion model with 10% CO ice to CO2 ice conversion. The chemical model includes a long prestellar core stage and a FHSC 
stage. The time axis is compressed when t < since evolution is slow during that time. 
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Fig. 10. — Top panel: The luminosity evolution of the episodic accretion model. The dash-dot line is the luminosity of the continuous 
accretion model for comparsion. Second panel: The total CO2 ice column density evolution in the episodic accretion model with 10% CO 
ice to CO2 ice conversion. Third panel: The gas phase CO column density evolution in the chemical model including the episodic accretion 
model with 10% CO ice to CO2 ice conversion. 
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Fig. 11. — The column densities of CO2 ice (total) from observations and each of the four luminosity-chemical model combinations. The 
upper panels have continuous accretion while the lower panels have episodic accretion. The left panels have standard chemical models 
while the right panels include the extra pathway for CO2 ice formation. The observed column densities are plotted as green dots, while 
the model predictions are plotted as black diamonds. 
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Fig. 12. — The gaseous C ls O column density from each of the four luminosity-chemical models is plotted versus time. The upper panels 
have continuous accretion while the lower panels have episodic accretion. The left panels have standard chemical models while the right 
panels include the extra pathway for CO2 ice formation. The dashed line is the typical observed column density of C ls O gas. 
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Fig. 13. — The gaseous CO abundance profile at 60000 yr from each of the four luminosity-chemical models. The time step is the last 
time step of the continuous accretion model. At this time, the episodic accretion model is in a high luminosity state, and the region of high 
gas phase CO abundance extends to much larger radii than in the continuous accretion model. However, the gas phase CO is still depleted 
in the inner regions due to conversion into CO2 ice. 
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Fig. 14. — The column densities of C ls O gas from observations and each of the four luminosity-chemical model combinations. The upper 
panels have continuous accretion while the lower panels have episodic accretion. The left panels have standard chemical models while the 
right panels include the extra pathway for CO2 ice formation. The observed column densities are plotted as green dots, while the model 
predictions are plotted as black diamonds. 



